These are fancy BEGIN plotmaps for taxonomic harmonisation. They are the data 2007-2010 and 2012-2014. 2011 is missing.
for (sp in sort(unique(BeginThinDraft$Species))) {
h <- BeginThinDraft %>% filter(Species == sp) %>%
mutate(Cover = as.numeric(Cover),
Plot = as.numeric(Plot),
Block = factor(Block, levels = rev(LETTERS[1:5])),
Cover = if_else(Cover == 0, NA_real_, Cover)) %>%
ggplot(aes(x = (Plot - 1) %% 3, y = (Plot - 1) %/% 3, fill = Cover)) +
geom_raster() +
facet_grid(Block ~ Year) +
scale_fill_viridis_c(limits = c(0, 100)) +
scale_x_continuous(expand = c(0, 0)) +
scale_y_continuous(expand = c(0, 0)) +
theme_classic() +
theme(axis.title = element_blank()) +
ggtitle(sp)
#ggsave(paste0("Maps/",i,".png"))
print(h)
}